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RADIO ASTRONOMICAL SPACE SYSTEM OF APERTURE 
SYNTHESIS! PILLING OP THE SPATIAL FREQUENCY 
SPECTRUM 

By 

N. S. Kardaishav^ S. V. Pograbanko, G. S.Tsaravskiy 

Introdu ction 

Recently more and raora attantion has bean given to the development of /3* 

long-range space radioastronomical ssystams (see [l,2j and references there, 

as well as C3»^3)* Space radioastronomical systems afford the possibility of 
practically unlimited increasa in the coilacting fstaas of the antennas in com- 
bination with the possibility of synthesizing high frequency two-dimensional, 
and even three-dimensional (]l3 radio Jjnages with resolution that is inaccessible 
under ground conditions. 

In this respect, it is important to examine certain informational poten- 
tialities of the radio astronomical space system of aperture synthesis Abbrevi- 
ation RASSAS) consisting of ground-based (GRT) and space (SRT) radio telescopes, 

i* Statement of Task 

One can demonstrate that the potentialities for the synthesis of radio 
images with RASSAS type systems are restricted by the sensitivity and resolution. 

As is known, the flux density from the source that has a brightness temperature 
of T and solid angle Q. is 

D 

F • 
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where k — Boltzmajin constant, K — wavelength. The root-mean-square value of the /k 
maximum detected flux from a point source during ohservations on a two-antenna 
radio interferometer is 

where T-- temperature of system noises, A — effective antenna area (A^rr D^/8, D« 
(D^xD^)^, D, and D — diameters of interferometer antennas), ^v— width of band of 
received frequencies, 6 t— integration time. After taking as the solid angle 
Tr(X/B)^/4 of the resolution of an interferometer with base B, and considering that 
P^^we obtain 



In the case where the antenna temi)orature, governed by the source is greater 
than the noise temperature of the receivers (strong source;, the inequality ( 3 ) 
becomes ^ 

^ ’ < 4 > 

where 0, F — .Angular size and flux of entire source, 6g, Pg — angular size and 
flux of component in the synthesized image. 

Assuming i6v®'2 MHz, At=lo3 s, Tg“10^^ K (limit for synchrotron omission), T= 

10 K, D*®30 in, 0=100 ang, s, P=100 Yan (1 Yan=l0“26w m"^ Hz”^), wo find that both 
for the weak and for the strong sources the limit base is about 2 million km. This 
indicates the outlook for the development of the RASSAS type systems that have 
the iiposslbillty of obtaining an extremely high angular resolution. 

It is important to examine the possible beginning variant for the space radio 
astronomical system (wo will call it RASSAS-l) that includes an SET that orbits 
around the earth on a low (H^=300-400 km) orbit. In light of the high relative 
rate of movement of the ground-based and space antennas of such a radio' .interfero- 
meter, a comparatively rapid filling of the spectrum of spatial frequencies is 
possible, and the diameter of the equivalent synthosizablo aperture is on the order 
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of the earth's diameter. From the ratio between the maximim size of the base B, 

diameters of the antennas GET and SET D and D , and the minimum detectable root- 

g s 

mean-square value of the brightness temperature in the image element T^i 
we obtain for the values: T=100 K, B=10^ km, D “30 m, D “100 m, 6v=2 MHz, At« 

n o & 

1 day— the amount T=10' K, This means that the synthesis system with low- 

a 

orbital SET can bo successfully employed to synthesize the images of compact sources 
of the typo nuclei of galaxies, quasars, stars, etc. Practically all the known 
quasars and galactic nuclei (total number of several ndrods) have that sur- 
pass the indicated limit, while unarbiguous radio images with resolution on the 
order of lO"^ ang, s, realized on interferometers with bases on the order of the 
earth's diameter, have boon obtained only foil several similar objects. 

In the examined system, the ahtennas participate in a complex relative move- 
ment that generally speaking is aporiodical, which makes it possible to obtain 
a fairly complete spectrum of spatial frequonciesj this significantly affects the 
unambiguity and quality of the obtained Images [5,, ll]. A similar' problem, in par- 
ticular, is encountered in processing and interpreting the data obtained by the 
ground-based in-tor for ometers that are used for synthesis of the earth's rotation, /6 
The latter are characterized by the elliptical shape of the motion tracks of 
the base vector projection in the plane of spatial frequencies (UV-plane) with a 
period of 12 hours [6]]. Since during each subsequent period the track is repeated, 
then the porcontago of filling of the UV-plano by the ground-based intoforomotor 
that consists of a small number of antennas is very small (see below), In con- 
trast to this, in EASSAS-1 practically complete and single-binding filling of 
the correspondii'»g region of the UV-plane occurs j its boundaries are close to an 
ellipse . 

In order to make a detailed study of the process of filling of the spatial 
frequency spedtrum in the RASSAS-1 variant, and to compare the results of filling 
with the filling done by the ground-based systems of aperture synthesis consisting 
of two (type VIB) and many (type VIBA) antennas, a series of programs were compiled 
for numerical modeling (in Fortran IV language, realization on ES-1040 comput;-^)* 
Below is a brief description of the algorithms (section 2), results of the 
modeling and their discussion (section 3)« A description and listings of the 
programs are given in the appendix. 
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2. Poser Ipt Io n of ?lc>t? >allng ithms _ 

2,1,' Coordinate Systg /m and Motion E quati ons 

The mevomant' of components In tha system of aperture synthesis Is examined 

in tha following rectangular coordinate system OXYZ that is generally accepted both 

$ 

in tha daecriptlon of radio intorferometars, as in the description of movement of 
artificial earth Batallites: tha center of the coordinate system "0** is located in 
the canter of the earth, tha axis ”0U" is directed towards the North Pole, the 
axis "OX'’’’ — towards tha point of spring equinox, and tha axis "OY" supplements the 
system of coordinates to the laft sat of throe, 

Tha motion aquations of tM» ground-based radio telescopes in the adopted system /7 
of coordinates look like; 

where R — earth's radius, XI —angular velocity of earth's rotaliion, ^ — geo- 

0 

graphical latitude and longitude of i-th GRT, — ajigle between the point of tha 

spring equinox and tha zero meridian at the beginning of the time reading. 


The motion aquS,tlons of the SRT on the assumption of a circular orbit and in 
tha absence of precession and nutation of the orbit look like 

, 6‘KH)*l^[C9nwt+e,)CoiSj, -Cott ^ 


(7) 


where R — ^radius of the SRT orbit, w— angular velocity of the SRT movement in tha 

0 

orbit, 0^ — initial phase of movement of SRT, XI ^ — longitude of ascending angle of 
orbit, i — wangle of incline of orbited plane to plane of earth's equator. 


Consideration for the affect of pi'oeassion of the SRT orbital plana on the 

result of filling tha UV-plana can be made by introducing tha dependence of o on 

time although, as will be shown later, for duration of observation T 1 day one 

o " 

can not consider tha precession. 


The motion aquations for the RASSAS-1 base looks like 


while the motion equations of the base (i, j)-the pair of radio telescopes of 
the ground-based system of synthesis look like 



(9) 


Sample geometry of the mutual arrangement and movement of the components in /8 
the RASSAS-1 system is presented in figure la. 


2.2. Observ ed Obje ct and Basis on UV -Plan e 


The position of the observed object in the celestial sphere is assigned by 
the coordinates a (direct ascent) and 6 (inclination). The direction vector to 
the object A^ is defined as 


a* 


‘Coifli CeiS, 

‘ Cm9, 


tlO) 


The observation plana (spatial analog of the UV -plane) is assigned, like the 
plana passing through the center ”0” of the coordinate system and perpendicular 
to the vector A. The basis vectors on the observation plane, or what is the same 
with accuracy to scale, on the UV -plana, are selected as equal 

•Coifli'CoiCSt’Vi) t 
-Colts' * 1 ^%) , 

Bvfc* &iv< CS ) , 

fcu ■’ C Bv ] • 

1. e.j the vector B is formed fxota A by turning it 90® in* the plane passing 
through A and the axis "OZ." The vectors A, Bj,,B.tfform an orthonormal basisj here 
the projection of the base vector in the system on B ^ determines the resolution for 
the coordinate 6 in the celestial sphace, the projection on B, —for the coordinate 
Ot, while the projection for A is the difference in the signal course from the source 
to one of the radio telescopes in relation to the other. The plan for selecting 
the basis on the UV-plane is given in fig. lb. 



ai;> 
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2.3. Calculation of SimuH 0^ 

For the CRT the condition of its visibility of the object corresponds to the 
fact that the radio telescope is located from tiie same side of the observation 
plane as the observed object, and for this the following condition must be fulfllledt 

.J3, 

For the SRT this condition must be supplemented, since depending on the 
height of its orbit the SRT will "see" the object in the space of a certain time 
after Its switching to the side of the observation plane that is opposite in rela- 
tion to the object until it enters the cylinder of shade created by the earth. 

Thus, the condition for cessation of the sparse radio telescope's visibility of the 
object can be written in the form 

In the process of simulation the filling’ of a UV-plane by a pair of radio 
telescopes is interrupted if even one of them at the given moment is not in the 
condition to see the object. 


2.4. Filling of the UV -plane 

In the course of the modeling the values of the projections of the current 
base vector of the given pair of radio telescopes on the baisls vectors B^and 
are determined. The obtained values are multiplied by the scale multiplier, 
selected such that the maximum base realizable by the system is written in the 
best manner into the discrete grid with assigned number of cells set aside for 
this. The projection values are rounded off to a whole number according to the /lO 

adopted quantization spacing. From the obtained whole-number projection values 
for the base vector on^,, and B^^ the nurabeys of the cells are determined in the 
computer storage files that have been set Jiside for storage of the results. 

The model result of a single observation at the given base during the time cor- 
responding to the adopted quantization spacing in time (in the given case, 
directly the quantization interval duration) is added to the contents of a cell 
thus determined, that is also symmetrical to it in relation to the center of the 
coordinates on the UV-plane. For a multiple -component system, this procedure is 
carried out in a cycle for all pairs of radio telescopes in the system. Thus, in 
the cells of the discrete presentation of the UV-plane the times for the stay of 
the system base projection in the corresponding areas of the continuous UV-plane 
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The deecrihod procedure also provides for the possibility of obtaining charts 
for the filling of the UV-plane during multiple-freciuency operation of the radio 
intetfaroroeter , The essence of the multiple-frequency filling consists of the 
fact that at the same time that a cell is filled that ha$ coordinates (u^, Vj) on 
the UV-plane corresponding, say, to the first ahd the highest frequency of obser- 
vation, colls ate also filled with coordinates ^ where < I 

is the -t-th relative observation froquancy. With the optimal selection of 
relative frequencies the filling of the UV-plane n,\ust be significant improved, 
of course, on the assumption that the distribution of radio brightness for the 
source remains unchanged in the entire range of observation frequencies. The UV- 
plane that is thus filled represent sbhe sum of the filling in the first frequency^ 
with the fillings obtained from it by the affine compression with the corresponding 
coefficients 'fl^ . 

Besides the increase in the number of filled calls on the UV-plane, the 
multiple-frequency operating pattern of the radio interfaroroetar has an important /ll 
qualitative advantage over the traditional one-frequency pattern that consists of 
the possible determination of the phase of the visibility function from the dif- 
ferences in the interference signal phases for each frequency pair. We propose 
to examine this question in detail in another- work, 

2. Processing of the Modeling Results 

The filling data obtained during the numerical modeling are further pcocessedj 

the following are determinedj value of the number of filled cells N^, maximum 

base realized in the given; numerical experiment, coefficient of time use 

(l.a., the relative portion of time during which simultaneous visibility of the 

studied object by both radio telescopes is realized), K.. In addition, the 

x 

function 'iJCt ) that is normed for a whole number of filled cells is de-terminedi 
c 

this is the function for the distribution of the number of filled cells according 

to the time of accumulation of t on one coll. 

c 

In the analysis of the different interferometric systems, it is convenient 
to select their information content as the comparison parameter. The information 
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content Sh according to Shannon is defined as 

Sk » S , <1;» 

w 

where (P/o^)j| j ratio signal/noisa in the (i, j)-th call of the discrete 

UV-plana, /i?;!, ^^»r--linaax dimensions of the call. With P/cL<l (weak source) one 
can be limited to the first term of the logarithm expansion into a Taylor series, 
and, taking into consideration that we obtain 

(K'.W)' ) . 

In the case of a strong source, p/c^» 1, the logarithm in (I5) can be /l2 
ignored as compared to the linear dapandanca of the information content on the 
number of filled calls, Thus, two estimates are posflble for the relative infor- 
mation content of the interior omatrlc systomsi according to formula (I6) for weak 
sources, arid according to the number of filled colls for strong ones. 

In this respect, during the processing of data from filling the UV-plana, wo 
determined the coefficient of information content Sh^ as well, 

2.6, Block Diagram of Counting and Presentatio n of Results 

The general block dia^am of the procedaro for numerical modeling of the 
process for filling the UV-plana is given below, 

DurintJ the printing of the program the initial data are issued and the 
results of processing the data on filling of the UV-plano, The charts for filling /I3 
of the UV-plana are derived on the ADP (alphabet-digital printer) in two gradationsi 
unfilled cells are represented by problems, filled— by the same symbol regardless 
of the time of accumulation on the given call, 

3* Results of Numerical Modeling and Their Discussio n 

3.1. Selection Initial Parameters 

A system was selected as the modal for RASSAS-1 that consists of one GRT 
that has the geographical coordinates of 

latitude 450 
longitude X«0°, 
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| Begin nlne| 




Input of Init ial data and constant b i 




s datarmlnatlon of tasis on UV-plantk! 


JkiU . 


baginning of cycla to model systam movaroanti 


J/ 


->“!datarmlnatlon of cursjant coordinates;— - 


V- 


'^iVarlf icatlon at entranca to shadow 


■ ■ 1. 


filling of currant call on UV-plana i 




and of cycla of modallng 


I 


'procassing of results i 


•output of modeling results and prlntingj 
r a nd , 


and one SRT that rotates around tlK,> earth on a circular orbit with tha following 
paramatarst 

altitude of orbit 

pariod of revolution 92 min, 

incline of orbital plana 
longitude of ascending loop 
baginning phase 


H «350 km, 
o 

T * 
o 

i= 52°, 


A «o®, 

o 

0 »QO 
o 

The geometry of tha system is presantad in figure la. 


It is evident that if tha duration of the observation pariod T^>T , then tna 

dapendanca of tha resulting filling of the UV -plane on 0 will be weak. With 

o 

24 hours the dependence on X will also be weak. For this reason, the zero 

values of these parameters ware selected purely formally. As will be shown below, 

the dependence o” is manifest as tha dependence on the difference of 

where a is tha direct ascent of tha observed object. In this respect, variations /l4 

were not made in tha value A . Tha values i and H ware selected close to tha 

o o 

corresponding values of the orbit of tha "Salyut” station. During tha numerical 
modeling calculations ware made for the situations that are distinguished by the 





FiguxQ ii Geometry of Radio Astronomical Space System of Aperture Synthesis 
(RASSAS-1) 

a) mutual arrangement of the ground-hased (GRT) and space (SRT) radio 
telescopes (B—earth's equator, G — Greenwich meridian) j b) plan for 
selecting basis on W-plane 

position of the observed object in the ctlestial sphere and the duration of the 
observation period, as well as by the number of analyzable frequencies during 
multiple-frequency filling of the UV-plane, Here, practically always, with the 
exception of the case whore the spatial response of the system was computed, the 
UV -plane was presented in the form of a mass of 101 x 101 colls f here the cell 
with the number (51,51) corresponds to the center of tljo coordinates on the UV- 
plano, i.o,, in our illustrations for clearness both halves of the UV -plane axe 
presented. The linear size AH of one cell is selected from the correlation 
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50 ^ 


ivn 

where R^--radlus of the SHT orbit, which corresponds to the fact that the 
maximum base realisable by the system aq,ualling will have a magnitude of 50 
cells in the discrato pEesantatlon. The spacing A t of the quantization of the 
process in time is sdected as equal to 6 seconds, which corresponds roughly to 
l/4 + l/lO of the average time that the projection of the system base vector 
stays in one call of the W-plane during a single passage through It. 

3.2. Dynamics of the Proces s of UV-D lane Filling 

Figures 2 and 3 present maps of filling for the UV -plane In the RASSAS-1 
system that ware computed without consideration for the precession of the SRT 
orbit for the observed objects with coordinates a“0^, 6*>45® (fig. 2) and (X“18*^, 

6*45® (fig. 3) with different durations of the observation period T^^. The numbers /l5 
in the left upper angles of the charts presented here correspond to the complete 
duration of the observation as follows 

No~"”~~TY' 12 13^4 5 6 7 8 

l^(hour)‘ i.5 I 3 j ^ I 12 ’ 24 ' 48 96 192 

On the charts corresponding to the short (T < 12 hours) periods of obsarva- 
tlon, the interruptions in the movement tracks of the projection for the system 
base vector on the UV-plane are clearly seeni they are due to the entrance of the 
SHT into the earth's shadow and the movement of the CRT during the time that the 
SRT is in the shadow, as well as the general epicycloid nature of the filling. 

Figure 4 presents graphs for the dependence of the number of filled calls 
(in thousands) on the duration of the filling period T^. The solid line repre- 
sents the dependence for the source with coordinates a**0^, 6*45*^^, and the dotted 
line- .for the source with a“18*^, 6»45®. It is apparent that, with the given 
quantization of the UV-plane the filling occurs fairly intensively aJl the way 
to time Tj^« 4 days, and then becomes saturated, which is induced by the repeated 
entrance of the projection of the system base vector into the previously filled 
cells. Hare the coefficient of filling a single-communications region that contains 
all the filled cells is fairly high, and is 0,5-0. 8. It is natural that the 
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Figura 2. Charts for Filling the UV-plana of tha RASSAS-1 Systam with: 

45 °, X“0®, H -350 km, T “92 min, i“52®,0 “0®, 0'-=O®, tt“ 

0‘\ 6“'!'5° for th8 filling t?ma T : 1.5 h (l)i 3 h (2), 6 h ( 3 ), 

12 h (4), 24 h (5), 48 h ( 6 ), 96 n (?), 192 h (6). Tha UV-plana 
is dividad into 101x101 calls. 

raduction in tha linaax dimansions of tha call, or what is tha sama, tha incraasa 
in thair number in tha discrata prasantation of tha UV-plana will raduca tha 
coafficiant of filling, and will shift tha momant that tha filling procass ba- 
comas saturatad toward the graatar timas. 

Pracassion significantly altars tha dynamics of UV-plana filling, but only 
at long obsarvation durations. Thus, with rata of pracassion 1-2 dag/day, tha 
linaar changas in tha basa vector with charactaristic siza R^“6700 km will ba 
200-400 km in 2 days tima, which corrasponds to 1-2 calls with the given 
quantization, Pracassion will raduca tha degree of covering of tha newly filled 
calls in the peripharal region of tha UV-plana, and correspondingly will in- 
crease their numi^u’, and, what is important from tha viewpoint of tha tcisks to ba 
solved by tha systam, on large values of the base. With T ^1 day, tha affect of 
pracassion cannot hava a significant effect on tha filling. 

The graph prasantad in figura 5 (solid line and dotted line correspond to 

the same positions of the observed object as in fig. 4") for tha dependence of 

the average tima of accumulation on ona cell t on the duration of the obsarvation 

av 

period T^ also indicates the effoct of repaatad entrance of tha basa vector 
projection into previously filled calls, and namely: up to time T^«2 days t^^ 
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Figure 3. Chart of UV-PlaJie Filling by the System RASSAS-1 with: 

.'i'=45o, X=0°, H ^=350 km, T^=92 min, i=52®,A^“0®, 6 q“0®, 

a=ie^, 6“45® for the filling time T : 1.5 h (1), 3 h (2), 

6 h (3), 12 h (4), 2'i h (5), 48 h (5), 96 h (?), 192 h (8). 

The UV -plane is divided into 101x101 cells, 

practically is not altered, which indicates the low percentage of repeated 

entrances, and then begins to increase rapidly (practically linearly with T„), 

n 

The effect of precession of the SRT orbital plane must slow this growth down at 
long observation times. 

3.3. Dependence of UV -Plane Filling on_ Coordinate of the Observe d Ob.ject 

The data presented in figures 2-5 indicate the significant change in the 

results of filling depending rn the coordinates (in the given case — on direct 

ascent) of the studied object. Thus, for the source with a*>18^, 6=45®, the number 

of cells filled in 24 h with other conditions equal, is roughly two times greater 

than for the source with a=0*^ and the same incline. This is linked to the fact 

that with different coordinates, the source"sees" the plane of the SRT orbit at 

different angles. The optimal from the viewpoint of the completeness of filling 

the UV -plane will be that mutual arrangement of the studied object and the SRT 

** ^ 

orbital plane, where the vector A of direction to the object will be perpendicular /I 7 
to the orbital plane} in other words, when the orbital plane will coincide with 
the observation plane. If one seaxcnes for the extreme only for a with the 
remaining parameters fixed, then it will be reached with -tt/2 for i<'rr/2 

and jQ^+tt/ 2 for tt/ 2< i < rr, where is expressed in radians. Due to the 

influence of both the SRT and CRT entering the shadow, and depending on and X, 
the value that* gives the maximum number of filled cells on the UV -plane for 
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duration of observation Tj^ (hours) 

Figure 4, Dejwndence of Number (in thousands) of Pilled Cells on UV-Plane 
on Duration of Filling^T„ withi 

^ X=0®, H =350 knu Tq= 92 min, i=52°, CL=0°, 0o"O«, 6-45°, 

a=18" (dotted line) a=0“ (solid line) with division of the UV -plane 
into 101 X 101 cells, 

the orbit with fixed pararae-* ors can be somewhat altered in relation to the value 
determined above. 

Figure 6 presents the charts for filling the UV -plane in 24 h for the objects 
with a=0*^ and 6 linked to the number of the chart as follows: 


No. 

1 

2 


4 

5 


7 i 

6 1 

0® 

i 

26® 1 

39° 1 

1^51° i 

64® 

77° 


8 

90® 


and figure 7 — the charts for filling the UV -plane in 24 h for objects with 6=45® 
and a linked to the number of the chart as follows; 


No 

1 

2 ! 3 i 4 

i i 

5 1 6 

7 

a ^ 

!o*^ 

J 

12 ^ j 15 ^ 

18^ 


8 


51 h 


The graph presented in figure 8a shows the dependence of the number of filled 
cells on the UV -plane (in thousands) on the incline of the studied object with 
cc»0*^ (solid line) and 0=18^ (dotted line), Tj^=24 h, and figure 8 b— the dependence 
of the number of filled cells on a of the studied object with 6=45® and other 

conditions equal. On figure 8b, the maximum with ce=18^ is clearly visible, which 
precisely corresponds to the case where the source with fixed 6 **seas” the plane 
of the orbit at the maximum angle. 

This fact has great importance for selecting the sequence of observation of 
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duration of observation Tj^(hoia:s) 


Figure 5, Dependence of Average Time of Accumulation on One Cell tgy on the 
Duration of Pilling Tj^ for* 

^=45°, H =350 km, T =92 min, i=52°, A^=0«, © =0°, 6=45°, a= 

0^ (solid line) and a=18*^ ^dotted line) Kith division of the UV-plane 
into 101 X 101 cells. 


several objects, since the presence of precession of the orbital plane in principle 
can permit the selection for each object of that moment of the beginning of obser- 
vations t , whore the filling of the UV-plano will bo the maximum possible' for^ it, 
o 

which will permit an increase in the efficiency of the observation process. 

The optimal moment for th-e beginning of the observations of the assigned object 
can be defined from the .correlation 

&oCto)“V4.=(X, ciB) 

The value t^ can be determined with accuracy up' to several days, for during such 
a period the effect of orbital precession on the process of filling the UV- 
plane Is insignificant. 

If one observes the sourc not only with the optimal value of the direct 

ascent, but also with the optimal value of incline 6 , then besides the ma xi m um 

o 

filling of the spectrum of spatial frequencies, the minimum frequency of inter- 
ference of the signals will occur. In fact, in this case the projection of the 
rate of SET movement on the vector A of direction to the object will equal zero, 
and consequently the change in the geometric delay (from here the interference 
frequency follows) will be determined only by the rate of movement of the ground- 
based telescope. 

* 
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Figure 6, Kaps of Filling the UV -Plane by the System RASSAS-1 with: 

T =24 h, .1=450, \..Qo, H =350 km, T “92 min, i“52®. P. “0°. 
r“0o, 6=450 and differSnt a: 0° ( 1 ) , 13° (2), 26° ( 3 ), 39° ( 4 ), 
51°C5), 64°(6), T7°(7), 90°l8). The UV -plane is divided 
into 101 X 101 cells. 


In practice a certain range of angles for the incline will exist, in which the 
projection of the SRT velocity on the vector A' will be fairly small, for example, 
it will not exceed the relative rate of movement of the radio telescopes in the 
VLB I type system with basis equal to the earth's diameter. We will introduce 
into the examination the amount of deviation 6 of incline from the optimal value: /19 
^6= 6-6 , Then the maximum value of the velocity projection of the Sf^T on the 
vector A can be estimated by the amount 

II I o 


Correspondingly, the maximum value for the projection of the relative rate of 
movement of the radio telescopes in the VLUI system will be 

U^,(VLfct) ^ ^ » ('^0. 

where T “24 h. Thus, one can estimate the range of angles at which V__(SRT)< 
Q II 

V^j (VLBl) from the correlation 
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which yields 
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or 


• /“ fc ii't' '* 

That Is, with the optimal a in the range of inclinations (6^-7®, 6 ^+ 70 ) the frequency 
of interference for RASSAS-1 will not he greater than during observations on the 
VLBI with base close to the earth's diameter. Consequently, in this case the 
requirements for the fast-response and the range of indefiniteness selection for 
the Doppler shift in frequency during observations on RASSAS-1 will be the same 
as for the ground-based system. In addition, an important consequence of selecting 
the optimal mutual orientation of the studied source and the orbital plane of the 
SRT is that the errors in determining the coordinates and velocity of the SRT, 
that are mainly linked to indefiniteness of its position along the orbit enter the 
analysis of the weight factor sin Ad, i,a., their influence will be significantly 
reduced , 

The value of the pure time synthesis is important from the viewpoint of /20 

efficient observations, i,e,, that part of the complete observation time during 
which simultaneous visibility of the studied source by both radio telescopes in 
the system was realized. Figure 9 presents the dependence of the pure time 

synthesis T on the coordinate a of the studied object with 6"45® and complete 

c 

observation time Tj^=24 h. It is apparent that with the value a“l8*^ that is . 
the optimal from the viewpoint of the complete filling of the UV-plane, the 
maximum pure synthesis time is also realized, 

3 , 4 , Distribution of t he Number of Filled Cells of the UV-plane According to the 
Accumulation Time 

With the adopted algorithm for filling the cells of the UV-plane, the 
accumulation times for one cell are distinguished among themselves, which is 

linked, first of all, to the: nonuniformity in the rate of movement of the projection 
for the system base vector over the UV-plane and secondly, to the repeated 
entrances of the base vector projection into the already filled cells, which occurs 
most often with small values of the current base. The range of change in time of 
accumulation for one cell is 2-'3 orders for one realization of filling? however, 
such' is. broad range is mainly due to the presence of cells with very large 
accumulation time? as a rule there are few of these (units). 
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Figure ?. Charts for Filling UV-Plane of RASSAS-1 System with: 

h, ^“45®, X=0®, H^=350 km, T^=92 min, i=52®,n^=0®, 

00=0®, 6=45® and different a; 0*^ (l), 3*^(2), 6*^(3), 9^(4), 

12^ (5), 15^ (6 ), 18^(7), 21*^ (8). The UV -plana is divided 
into 101 X 101 calls. 

The graphs presented in fig. 10a, b show the functions ([) (to) for the dis- 
tribution of the number of filled calls according to accumulation time on one 
call t , obtained for T ®24 h, with coordinates of the studied object a=18*^, 6^5®» 
and division of the UV-plane into 101 x 101 (fig. 10a) and 128 x 128 (fig. 10b) 

cells. The maximum values of t plotted on the horizontal axes of the given /Z1 

c 

graphs correspond to the maucimum values for the accumulation time for one cell 
that ware realized during the given numerical experiments. Both graphs have 
clearly pronounced maiximums at accumulation times close to the average time of 
accumulation t^^ shown on the graphs by a vertical line. It is apparent that 
with an increase in the number of colls in the discrete presentation of the UV-plane 
the maximum of the distribution function appears more distinctly. The percentage 
of points with accumulation time lower than 0.5 t in both cases is small («10%), 
and is reduced with an increase in the number of colls in the division of the UV- 
plane, 

It is worth noting that the cells with large accumulation time can be used 
for the purpose of calibrating the system and for accurate phase correlation 
in solving the tasks of radio astrometry. 
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Inclination 6(°) 


a 



dizact ascent a (n) 


b 

Plguro 8. Dapendance of Nuxnbar of Fillad Galls N„ (in thousands) on UV- 
Flana on Inclination 6 of Studied Object for* 
a=0" (solid lino), and a»»18" (dotted lino) with 24 h, 

i = 450 , \=0°', H =>350 km, T =92 min, i«52°, D. -0°, e-O®. 

The UV-plano is°dividod in?o 101 x 101 colls,® 

b. Dopondonco of Number of Pilled Colls No (in thousands) 
on UV-plano on Direct Ascent of a of Studied Object with 
6*»45° and ^o other conditions the same as in fig, 8a, 



),5» Mult Ipla -Freque ncy Filling of UV -Plane 


Figure 11 presents charts for filling of the UV -plane In 24 h for the source 
with a"18^, 6 * 45 ® with one-, two-, three- and four -frequency filling with 
relative frequencies 1.0, 0.9» 0.8 and 0.7, and figure 12 shows the dependence 
of the number of filled cells on the number of filling frequencies for the 
sources with a“0^, 6“45® (solid line) and a*=18*^, 6*45° (dotted line). The values 
that are discrete in their essence, presented in fig, 12, are connected by a 
smooth curve for clarity. 

On figure 12'; one can see the significant growth In the number of filled 
cells during the transition from one frequency to two, and the saturation in the 
growth during further increase in the number of frequencies. However, with 
varying spacing of the difference in relative frequencies, different stages in the 
filling increase are possible. In particular, the charts presented in figure 11 
show the possible additional increase in the number of filled cells during the /22 
introduction of yet another frequency, closer to 1,0 than 0.9, since this would 
permit the filling of the gaps preiwnt in the peripheral regions of the examined 
charts I the more so since this gives the advantage to the number of cells that 
have high values of spatial frequencies, 

3,6, Spatial Response ^ the RASSAS-1 System 

The spatial response (beam pattern of the synthesized aperture) in the 
RASSAS-1 system was obtained by a Fourier transform of the UV-plane filling 
obtained during modeling of the observation of a source with a®18*^, 6*45°, Tj^** 

24 h, and with the division of the UV-plane into 128 x 128 cells. 

The obtained response, constructed in the form of a relief map, is presented 
in figure I 3 . The maximum lateral lobe of the given response is 12^ of the 

magnitude of the central maximum, while the root-mean-square level of the lateral 
lobes equals 1,2^, The maximum and the root-mean-square levels of the lateral 
lobes can be significantly reduced, both with an Increase in the duration of the 
observation T^^ above 24 h, and with the use of the extant methods of image 
clarification [8,9^, 



dlroct ascent a (n) 


Figure 9, 


Dependence of Pure Synthesis Time I 
Ascont a of Studied Object with* 
T„»2^ h, X-Oo, H «350 km, T » 

6845 ° ° ° 


in RASSAS-l System on Direct 
92 min, i=52®,A^“0®, 6*0®, 


3.7. Comparison with the Ground-Based Systems Type VI£I and VLBA 


To compare the filling of the UV-plane by the RASSAS-1 system with the filling 
by the created ground-based systems of aperture synthesis using the earth's rotation, 
two systems were selected of radio Interferometerst a) two-component VliBI with 
antennas that have the following geographical coordinatesi 

No of antenna latitude longitude 

1 42.0® -80,0® 


44 . 5 c 


34.0‘ 


/23 


which roughly corresponds to the interferometer Simeyz-Khaysteki b) the VLBA system 
suggested in [[loj, consisting of 10 radio telescopes whose geographical coordinates 
are given below 


of antenna 

Latitude 

Longitude 

I 

40.4® 

- 3.2® 

2 

38.4® 

- 77.2® 

3 

S3, 6® 

- 79.8® 

4 

38.3® 

- 81.5® 

5 

38.5® 

- 06.6® 
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Figure 10. Distribution of Number of Pilled Cells on UV -plane according 

to Accumulation '^line for One Call t vrltht 

T„“24 h, 'l?=45o, X=0°, H “350 km, T°“92 min, i“52®, A -0°, 
il L o o o 

0^®^°, 0C“18 , 6“45°, and with division of the UV -plana intoj 
a) 101 X 101 cells, b) 128 x 128 cells. 
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P'igura 11, 3harts for Filling tho UV -Plano by the RASSA2-1 Systom with Multlplo- 
Froquancy Filling with Rolall ''0 Fxaquanciasi 

-r) ■= 1.0, 0.9, 0.8 and 0.7 with T^»24 h, A.-0®, H "350 km, 

T^"92 min, i"52°,0^*^°. Qq-O®, a-18^^, 6“45®, division of tho UV- 

plano into 101 x 101 colls and with varying numbor of fraquanclos 
n“l, 2,3,^. 

It was consldarod in tho calculations that tho dimonslons of all tho antonnas 
in tho systems RAGSAG-1, VLBI and VLBA wore tho samo, Tho scalos of tho charts 
for tho UV -plana for tho VLBI and VI£A systams wara soloctad from tho corrolationsi 

60 ^ , 

50 a 2 ti# 

in ordar to raduco tho valuos for tho maximujii banos realizabla by this systam, 
axprassad in tha numbar of UV -plana calls, to tha valuos roallzabla by tha 
RASSAS-1 systam, Tho obtained valuos fiR axo 212 km for VLBI and 256 km for VLBA, 
which is closa to Ar«272 km for RASSAS-1, therefore tho charts foi filling of 
tho UV-plana and thair statistical characteristics obtained for different systams 
can bo directly compared among themsalvas by taking into consideration the fact /24 
that tha numbor of antanaas participating in tha VLBA operation is considerably 
greater than the, numbar of antennas in tho RASSAS-1 system and tho VLBI systam. 

Figure 14 pcvisants tha charts for filling of tho UV-plana by tha VI£I 
modal system for objects with 6 corresponding to tha chart numbor as follows* 


No 

1 

2 

3 

4 

! 5 

6 

6 

0 ° 

^15^ 

30® 

450 

60 ® 

90 ® 


Figure 15 presents charts for tha filling of tha UV-plana by tho VLBA 
systams (index "G" — ground-basad at tho chart numbar) and RASSAS-1 (index "S" — 
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numbar of f*raquenclas of filling 
UV -plana 

Plgura 12, Dapendanca of Numbar of Pillad Galls N« (in thousands) on UV -Plana 
with Multlpla-Fraquancy Pilling on Numbar of Simult^aously Ana- 
lyzabla f'raquancias with Ralatlva Praquanciast 
'^^"1.0, 0,9, 0,8, 0,7 and Tj^«24 h, X«0*, H =350 min, T.«» 

92 min, 1«52®,^"0®, 0 “O®, 6-45®, a»o" (solid lina) and 0 t“’ lo ^ 
(dottad lina) , o 

UV -plana dlvidad into 101 x 101 calls, 

spaca)for objacts with inclination corresponding to the chart numbar as followst 


No 

1 

2 1 3 j 

4 

6 

0 ° 

45® 1 60 ® j 

90 ® 


For the inter far omatric systems type VUBI and VLBA the value of the direct 
ascent of the studied source with daily synthesis is quite insignificant, and 
for the RASSAS-1 in the given numerical exparimants d was selected as equal to 
18. 


For comparison of the results of filling the UV -plana by the model systems 
VLBI, VLBA and RASSAS-1, tha table presents the values for the number of filled 
cells N^, the relative information content for a weak {Signal Sh^, and the maximum 
base B realizable in the given numerical experiment. The complete observation 
time in each case equals 2^'f h. The table corresponds to the situations whose 
charts are presented in fig, 14 a^id fig. I 5 . 

Conclusions 

The calculations presonted above, in particular, the charts for tha 
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filling of the UV’ plane and the data presented in table 1, make it possible to 
draw the following conclusions: 

1. The filling of the plane of spatial frequencies by the RASSAS-1 system 
that consists of two radio telescopes (one ground-based and one space) is much 
higher than the filling done by the ground-based two-component intercontinental 
radio interferometers of the type VLBI, 

2. Pilling of the UV-plane by the RASSAS-l system in Zk hours is signi- 
fj.cantly better (on the average) than the filling done by the multiple-antenna 
(up- to 10 antennas) global interferometric systernsj here with an increase in 
the dvwation of observations over 24 h RASSAS^'l makes it possible to obtain 
even more complete filling, which, of course, is impossible in the VLBA system. 

3. In the case of observance of a weak source for 24 h, the information 
content (according to Shannon) of RASSAS-1 on the average is three times inferior 
to the information content of VISA (10 telescopes), which is evidently due to the 
considerable differences in the collection areas. It is easy to show that the 
indicated Information content is compared for both systems with the 7-anterma 

26 





Figure 15. Charts for Filling UV-Plajie in T =24 h by VLBA Systems (index "G" 
and RASSAS-1 (index "S") for Studied Objects with 0=18^, 6=0® (l) 

45® (2), 60® (3), 90 ® (4). Division of the UV-plane into 101 x lOl 
cells 

VLBA variant. But such a system will fill already roughly three times fewer cells 
of the UV-plane as compared to RASSAS-1. 

4. In the comparison of the time-linked observation program for the RASSAS 
type systems it is necessary to consider the direct ascent of the sources and 
the precession of the SRT orbital piano in order to obtain the optimal filling 
of the spectrum of spatial frequencies for each source included in the program. 

5. The described algorithms and the programs of numerical modeling for the 
filling of the UV-plane by the RASSAS-1 system can be used +0 select objects for 
study and compilation of the time sequence of their observations in real 
experiments on a space radio interferometer, 

I^eferencas 

1. Buyakas, V. I.j Gvamichava, A. S.j et al. , Neog r an ichenno n ar ashc hiva yemy y 

kosmic h~»skiv radioteleskop L”Unlimited-Accumulation Space Telescope'^, Pre- 
print of the Institute of Space Research, USSR Academy of Sciences, No 373, 

1977. 

2. Buyakas, V. I.j Gvamichava, A. S. j et al. , Kosmic hesklye issledovaniya, 

16 (1978), No. 5. 

3 . Powell, N. V.{ Hibbs, A. R. As tron. Aeronaut. 15, No. 12 (1977). 58. 


27 


4. Easier, R. P.f Johnson, G. L.| and Vondrak, R. N, Radi o ScjU, 12, No. 5 
(1977). 845. 


5. Fomalont, B.| and Rayt, M, in G alakt Icheskaya 1 vnegalaktlche ska^a. 

r adloastr onomlya ["Galactic and Extragalactic Radio Astronomy"^ Moscow, 

Mir, 1976. 

6. Yasepkina, N._ A. j Korol'kov, D. V.| and Pariyskiy, Yu. N. Radiotelaskopy i 

radiometry ["Radio Telescopes and Radiomotry"3» Moscow, Nauka, 1973. 

7. Aksenov, Ye. P. Teorlva dvizhanly a iskusstvennykh sputniko v Zemli ["Theory 

of the Movement of Artificial Earth Satoliltes"[|, Moscow, Nauka, 1977. 

8. Hogbom, J. A. Astron. Astrophys. , Suppl,, 15 (1974), 417. 

9. Wornecka, S. J., Radio S ci.. 12, No. 5 (1977), 83 I. 

10. Kellarmann, K. I. VLB I Network Studies . III. An Interc ontinental Very Long 
B aseli ne Arrajj NRAG Report, 1977. 

11. Kogan, L. V. Sver zhdal * nyaya interferometriy a kompaktnykh rad ioisot chnikov 

["Ultra Long-Range Interferometry of Compact Radio Sources"3r Candidate 
dissertation. Institute of Space Research, USSR Academy of Sciences, 1972. 


Appendix. Prograia of Numerical Modeling of UV-Pla ne F illin g by Radio Astro- 
n omical Spac e Sjj^stom of Aperture Synthesi;^ 

The prograjn of numerical modeling of UV-plane filling by a radio astronomical /42 
space system of aperture synthesis consists of the main program (MAINPGM) and 
the subprograms ANGL, oPBII, MZBRO, SCALE, BASIS, VPROD, STATE and OUTUV. 

The MAINPGM implements the input and control opening of the constants and 
initial data, and gives the initial data in the necessary form. The initial data 
and the constants are introduced with the help of the operators DATA and appro- 
priation, The constants are PI*»3*l4..., radius of the earth RE, acceleration of 
the gravity force on the earth’s surface G, angular velocity of the earth’s rotation 
WE. The initial datat PHl(l), PHI (2)— geographical latitude' and longitude of 
the GRTi PSI (1), PSl(2), PSI<3), H— longitiide of ascending angle of SRT orbit, 
incline of orbital plane, beginning phase at moment of time TO and height of or bit j 
ALP(l), ALP (2) — direct ascent and inclination of observed sourcef TO, T, DT — 
beginning of time reading, duration of observation, and spacing of quantization 
of process in time. All the distances are expressed in km, the time in min, the 
angles (except ALP(l)) in degrees, AIP’(l)— in hours. During its work the MAINPGM 
generates the subprograms ANGL, ORBII, STATE and OUTUV, 

The subprogram ANGL translates the angles from degree measures to radians. 

The input parameterss RAD — number of radians par degree, ANG — identifier of angle 
file subject to transformation, N — length of file. The output values for the 
angles are contained in the ANG file. 

The subprogram ORBII strictly models the process of UV-plane filling by the 
RASSAS system. The input parameterss PI=3.14,.,, P— identifier of file, removed /43 
for the UV-plane area to be filled, PSI, PHI, ALP — files of angular parameters of 
the system, N — size of the file P, NO — coordinates for the center on the UV-plane, 

RS, RE — radii of orbit and earth, RM — inverse size of one cell on the UV-plane, 

WS, WE — angular velocities of SRT and earth rotation, T, D, T — duration of 
observation and spacing of quantization with respect to time. After fulfillment 
of the subprogram, all the parameters except P .remain unchanged, and the cells 
of the file P contain filling of the UV-plane, i.e, , the tima for stay of the 
projection for the system base vector in the given cell of the UV-plane, The 
generated subprograms* lEERO, SCALE, BASIS. 
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Subprogram MZERO brings to zero the two-dimensional file X of size N x N, 

The subprogram-function SCALE makes a scalar ‘product of two three-component 
vectors presented in the form of files A and B of length 3 cells. 

The subprogram BASIS determines the basis on the UV-plaxie, The input 
parameterst ALF — file 2 cells long containing the source coordinates. 

The output parameters: A. BU, BV — files 3 cells long corresponding to the vector, 

direction to the studied object, and basis vectors on the UV -plane. The generated 
subprogram is VPROD, 

The subprogram VPROD makes a vector product of two three-component vectors 
A and B, The result is arranged in the file G of length 3 cells. 

The subprogram STATB statistically processes the results of UV-plane filling. 

The input parametarsi F — file of calls of tha UV-plane, N— its size, NG— 
coordinate of the cantor on the UV-plana, DR — size of one coll in km, T — duration 
of observation. During the operation of tha subprogram the following amounts 
are deter mined and issued in print: N — number of filled colls, NO — number of colls /44 

in circular region with radius equal to tha maximum realized base, KF— coeffi- 
cient of filling of circular region, BM — maximum realized base, XM and YM — 
maocimum realized values for tha base projections on the axis U and V, TMIN, 

TMAX, TBV — minimuiri, maximum and average time of accumulation for one coll, TUSE — 

"pure” synthesis time, KUSE — coefficient of use of the time in the given numerical 
experiment, AINF — information content of the synthesized aperture for a weak 
signal. In addition, tha distribution function for the number of filled cells 
according to tha accumulation time for one cell is determined and issued in print. 

An example of the output of the statistical characteristics is given below. 

The subprogram OUTUV produces a display in prln-c of the UV-plane filling 
in two gradations. The input parameterst F— identifier of derived .two-dlmansional 
file, N — its size. 

The listings of the described subprograms are given below. 
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00 |/|| fO«TtAK 




data a,yl/p,3lA»4.lJs>lt-J/ 
pATA MI/41. 0,0,0/ 

DATA 0||/0.A,|)f.0,0.0/ 

data ai)/i4.o.iI.o/ 

100 TOAHAT|//f», 1 lO# 1 H*» . 

♦ / / 1 Oil ' . ) 


1*0 » « «" * I w 1 J ! ' : -S' y , 

♦ y iox f ♦ 0 T»» iPf . 1 1 1 K » 




♦ 1 *,» ^»^IO*»*Ba« 

ITO Jo*«iM/ 10 K,fFMHO ) 


T W"W . W 

,5!,to*P<ie-i, /«»*•*. 0 ) 

Bl*1 .O/RM 

rSISt 

58 lu: . 

58 Si |J8: St . 

OR Rt 100,D* 

OR NT 1T0 
RAD«^t/li0,0 , 

CALL AN0L(RAD . 1 , 2 > 
(ALL ANOlNaD.RII . 3 ) 

call STiTRCF ,N,RC , 0 R»T) 
OuTUV<P»N) 

Imo 


r;>;’»i/l« rORT*An iv V 


iUIROUTINI AN0L<RA0 » ano ,r> 

0 l**INIION ANC(R) 

DO ^ I • 1 I N 

ARf ( 1) «AN0( I> *RAO 

Rf TURR 

f RO 
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